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ABSTRACT 

We propose a simple way to estimate the parameter /3 ~ f2°-^/& from three-dimensional 
galaxy surveys, where Q, is the non-relativistic matter density parameter of the Universe and 
h is the bias between the galaxy distribution and the total matter distribution. Our method 
consists in measuring the relation between the cosmological velocity and gravity fields, and 
thus requires peculiar velocity measurements. The relation is measured directly in redshift 
space, so there is no need to reconstruct the density field in real space. In linear theory, the 
radial components of the gravity and velocity fields in redshift space are expected to be tightly 
correlated, with a slope given, in the distant observer approximation, by 



g ^ y/l + 6/3/5 + 3/377 
V 13 

We test extensively this relation using controlled numerical experiments based on a cosmolog- 
ical A^-body simulation. To perform the measurements, we propose a new and rather simple 
adaptive interpolation scheme to estimate the velocity and the gravity field on a grid. 

One of the most striking results is that nonlinear effects, including 'fingers of God', affect 
mainly the tails of the joint probability distribution function (PDF) of the velocity and gravity 
field: the 1-1.5(7 region around the maximum of the PDF is dominated by the linear theory 
regime, both in real and redshift space. This is understood explicitly by using the spherical 
collapse model as a proxy of nonlinear dynamics. 

Applications of the method to real galaxy catalogs are discussed, including a preliminary 
investigation on homogeneous (volume limited) "galaxy" samples extracted from the simula- 
tion with simple prescriptions based on halo and sub-structure identification, to quantify the 
effects of the bias between the galaxy distribution and the total matter distribution, as well as 
the effects of shot noise. 
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1 INTRODUCTION 

Analyzes of large-scale structure of the Universe provide estimates 
of cosmological parameters that are complementary to those from 
the cosmic microwave background measurements. In particular, 
comparing the large-scale distribution of galaxies to their peculiar 
velocities enables one to constrain the quantity (5 = OP'^ /h. Here, 
Q. is the cosmological non-relativistic matter density parameter and 
h is the linear bias of galaxies that are used to trace the underly- 
ing mass distribution. This is so because the peculiar velocity field, 
V, is induced gravitationally and is therefore tightly coupled to the 
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matter distribution. In the linear regime, this relationship takes the 
form (Peebles 1980) 



(1) 



where 5 denotes the mass density contrast and distances have been 
expressed in km • s~^. Under the assumption of linear bias 5 — 
b~^5^^\ where 5^^^ denotes the density contrast of galaxies, the 
amplitude of peculiar velocities depends linearly on (3. 

Density-velocity comparisons are done by extracting the den- 
sity field from full-sky redshift surveys (such as the PSCz, Saunders 
et al. 2000; or the 2MRS, Erdogdu et al. 2006), and comparing it 
to the observed velocity field from peculiar velocity surveys (such 
as the Mark III catalog, Willick et al. 1997; ENEAR, da Costa 
et al. 2000; and more recently SFI-l~l-, Masters et al. 2005). The 
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methods for doing this fall into two broad categories. One can use 
equation (T) to calculate the predicted velocity field from a redshift 
survey, and compare the result with the measured peculiar veloc- 
ity field; this is referred to as a velocity-velocity comparison (e.g., 
Kaiser et al. 1991; Willick & Strauss 1998). Alternatively, one can 
use the differential form of this equation, and calculate the diver- 
gence of the observed velocity field to compare directly with the 
density field from a redshift survey; this is called a density-density 
comparison (e.g., Dekel et al. 1993; Sigad et al. 1998). The advan- 
tage of density-density comparisons is that they are purely local, 
but they are significantly sensitive to shot noise because the diver- 
gence of the observed velocity field is estimated from the sparse ve- 
locity sample. The integral form of velocity-velocity comparisons 
make them much less sensitive to such a noise, but their non-local 
nature make them affected by the tides due to the presence of fluc- 
tuations outside the survey volume (see, e.g.. Kaiser & Stebbins 
1991). 

Still, the common problem with both types of these compar- 
isons is that while equation ([TJ involves the density field in real 
space, the observed density field is given solely in redshift space. 
One approach to tackle with this problem is to reconstruct the real- 
space density field from the redshift-space one. The transformation 
from the real space coordinate, r, to the redshift space coordinate, 
s, is 

s = r + i'(r)f, (2) 

where v{r) = v • f ; v(r) is the real-space velocity field, and ve- 
locities are measured relative to the rest frame of CMB. Therefore, 
one has to correct galaxy positions for their peculiar velocities. To 
do so, equation (TJ is used and a self-consistent solution for the 
density field in real space is usually obtained iteratively (Yahil et 
al. 1991). In the first iteration, to predict peculiar velocities accord- 
ing to equation (TJ, the real-space density field appearing on the 
r.h.s. of this equation is approximated by the redshift space-density 
field, and so on, until convergence. However, from equation (TJ it 
is obvious that the amplitude of peculiar velocities depends on /3, 
the parameter to be subsequently estimated. Therefore, to perform 
a density-velocity comparison self-consistently, one has to recon- 
struct the real space density field for a range of different values 
of /3. For example, Branchini et al. (1999) performed such recon- 
structions for 10 different values of j3 in the range 0.1-1.0. Note 
that instead of this traditional algorithm, more sophisticated meth- 
ods now rely on Euler-Lagrange action minimization (e.g., Pee- 
bles 1989; Shaya, Peebles & TuUy 1995; Nusser & Branchini 2000; 
Phelps 2002; Phelps et al. 2006) or resolution of optimal assigne- 
ment problems (e.g.. Croft & Gaztanaga 1997; Frisch et al. 2002; 
Mohayaee et a. 2003; Mohayaee & Tully 2005). 

Another approach, proposed by Nusser & Davis (1994; here- 
after ND), is to perform the comparison directly in redshift space. 
ND derived the density-velocity relation in redshift space in the lin- 
ear regime. Because velocity-velocity comparisons seem in prac- 
tice more robust than density-density ones, they aimed at trans- 
forming this relation to an integral form, i.e. at solving for the ve- 
locity as a functional of the redshift-space density. Due to the radial 
character of redshift-space distortions, it turned out to be possible 
only via a modal expansion of the density and velocity fields in 
spherical harmonics. Using this expansion, ND were able to put 
constraints on [3 (ND; Davis, Nusser & Willick 1996). However, 
the approach with the reconstruction of the real-space density field 
remained popular. Apparently, equation (TJ is appealing by its sim- 
plicity, both in terms of reconstruction of the velocity field, and of 
estimation of the parameter j3. 



This paper is devoted to finding an equivalent of equation (TJ 
that would hold for redshift-space quantities, but would share its 
simplicity. Specifically, let us define the scaled gravity: 

g(r).|^g:5(«)(r')^. (3) 

Under the assumption of linear bias, adopted here, g is proportional 
to gravitational acceleration, and can be directly measured from a 
3-D galaxy survey. Equation (TJ then implies 

g = /3"'v. (4) 

Now, let us assume that we measure the gravitational acceleration 
directly in redshift space: 

g.(s)^/5^^<^'(s')^, (5) 

where 5s*^ denotes the density contrast of galaxies in redshift 
space. We will follow ND in a natural definition of the redshift- 
space velocity field: 

Vs(s) = v(r) . (6) 

What is the relation between gs and Vs? 

Equation (4j holds strictly in the linear regime. Nevertheless, 
numerical simulations (Cieciel^g et al. 2003) have shown that it 
remains accurate to a few percent for fully nonlinear gravity and 
velocity fields. These results will be fully confirmed by the present 
work and can be explained by the fact that velocity and gravity 
fields are dominated by long- wavelength, linear modes. Therefore, 
in deriving the redshift-space counterpart of Equation (4j, we will 
apply linear theory. Unfortunately, there is no local deterministic 
relation between gs and Vg. However, as shown below, gs and Vs 
are strongly correlated, so the mean relation will be a useful quan- 
tity. Only radial components of velocities are directly measurable, 
so we will be interested in the relation between gs = gs • s and 
t;s = Vs ■ s. 

This paper is thus organized as follows. In §|2l we compute the 
properties of the joint probability distribution function (PDF) of the 
fields Qs and Vs in the framework of linear theory and distant ob- 
server limit. In the linear regime, this PDF is expected to be a Gaus- 
sian, entirely determined by its second order moments. The quan- 
tities of interest are {vsgs)/{vl), {gl) / {vsgs) and {gi)/{vi}: 
they all give an estimate of the expected ratio between g^ and t^s 
and the difference between them can be used to compute the scatter 
on the relation, that will be shown to be small. The validity of our 
assumptions, in particular the distant observer approximation, is ex- 
amined a posteriori in §[3] We also justify our choice of CMB rest 
frame redshifts, needed to avoid the so-called rocket effect (Kaiser 
1987; Kaiser & Lahav 1988). In §[4] linear theory is tested against 
numerical experiments, both in real and redshift space. To do that 
we use a dark matter cosmological A'^-body simulation with high 
resolution, allowing us to probe the highly nonlinear regime. We 
propose a new and simple algorithm to interpolate the velocity and 
gravity field on a grid from a distribution of particles. We address 
extensively a number of issues, such as the validity of the distant 
observer limit, edge effects, cosmic variance effects, effects arising 
from non-linear dynamics, in particular so called 'fingers of God' 
(FOG) and effects of dilution (numbers of tracers used to probe 
the velocity and gravity fields). We use spherical top-hat model to 
support our interpretations of the measurements (technical details 
are given in Appendix A). The effect of the bias is also examined 
briefly by extracting from the dark matter distribution two kinds of 
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subsamples, one where each "galaxy" is identified with a dark mat- 
ter halo, and the other one where each "galaxy" is identified with a 
dark matter substructure. Finally, the main results of this work are 
summarized in §[5] In this last section, we discuss observational is- 
sues such as discreteness effects, appropriate treatment of biasing, 
incompleteness, errors on peculiar velocity estimates, etc. These is- 
sues will be addressed in detail in a forthcoming work, where the 
method will be applied to real data. 



2 LINEAR THEORY PREDICTIONS 



Linear regression of Va on gs yields 

^= - {vi} • 



(7) 



The symbols (•) denote ensemble averaging. We can thus charac- 
terize the linear relation between gravity and velocity (gravity in 
terms of velocity) by its slope. 



(8) 



Alternatively, one can study the inverse relation (velocity in terms 
of gravity). Linear regression then yields Vb — ((^'si7s)/{gs))i7s, 
so the inverse slope is aj = («s5s)/((/s )• To describe the linear 
velocity-gravity relation, we have thus at our disposal two estima- 
tors of its slope, the forward slope, ap, and the reciprocal of the 
slope of the inverse relation. 



I ^ / \ ■ 
Due to the scatter, these two estimators are not equal: 



where 



(9) 



(10) 



(11) 



is the cross-correlation coefficient of the velocity and gravity fields; 
f^v ~ {Va) and (Tg — (gs)- The linear regression, although the best 
among all linear fits to a cloud of points, visually looks biased. For 
two correlated Gaussian variables, an unbiased slope of the isocon- 
tours (ellipsoids) of their joint PDF is given by the square root of 
the ratio of the variances of the two variables. Thus we have third 
estimator of the slope. 



au 



(12) 



This is easy to show that it predicts intermediate values between 
aj^ and a f ■ 

Second-order moments of the joint PDF for gravity and ve- 
locity, which appear in the above three estimators of its slope, are 
much easier to perform in Fourier space. We will adopt the distant 
observer limit (DOL), in which the Jacobian of the transformation 
from real to redshift space, equation simplifies to 



J(r) 



1 + 



dv 
dr 



(13) 



The Fourier transform of the redshift-space velocity field, using the 
mapping ([2}, is 



Vs(k) = / e* ''vs(s) d^s = 



ik'r ikv / \ t/ \ i3 

e e v(r) j(r) d r 



(14) 



where v(k) is the Fourier transform of the real-space velocity field. 
The linearized continuity equation yields v(k) — — /(ik/fc^)5(k), 
where / ~ OP'^ and 5(k) is the Fourier transform of the real-space 
density field, hence 



From equation ([5} we have 
g.(k) 



(15) 



(16) 



This equation looks similar to the preceding one, but here appears 
the Fourier transform of the redshift-space galaxy density field. 
Moreover, unlike the preceding equation, the above equation is ex- 
act. 

Radial components of the redshift-space velocity and gravity 
fieldsareSs(k) = -/ (i/i/fc) (5(k) and3s(k) = -(i^i/fc) (k), 
where = k ■ s/{ks). From conservation of numbers of galax- 
ies in real and redshift space, it is straightforward to write down 
an equation for the Fourier transform of the redshift-space galaxy 
density contrast. In the linear regime it reduces to 5s*'' (k) — 
(5(s'(k) + f^?~5{^^) (Kaiser 1987), or 



(17) 



Therefore, we obtain 
Os(k) 
and 



(18) 



S.(k) = -6(l + /3M^)^5(k). (19) 

The above pair of equations enables us to calculate the aver- 
ages appearing in equations l[8j, l|9), and il2t . Specifically, 

{vD = (2^)-y ^-^(^+^'ys fitf'Z (5(k)5(k'))d^fcd^fc'.(20) 

For a homogeneous and isotropic random process, ((5(k)<5(k')) = 
(27r)^5i3(k + k')P(fc), where So is Dirac's delta and P{k) is the 
power spectrum of the real-space density field. Performing the in- 
tegral over k' yields 



p{k) dk. 



(21) 



Similarly, 



bf[l + (3/5)/3] 

/o 



67r2 



P{k)dk. 



(22) 



Finally, 

/ 2^ b 



67r2 



, 6/3 3/3^ 



P{k) dk. 



This yields 



(vsgs) 1 + 3/3/5 



(23) 



(24) 
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I) 



au 



M) 



1 + 6/3/5 + 3/3 V? 
+ 3/3/5) ' 

^1 + 6/3/5 + 3/377 



(25) 



(26) 



Equations ll24t-(l26t provide a way to estimate /3 compar- 
ing gravitational accelerations to peculiar velocities of galaxies di- 
rectly in redshift space. All three estimators predict the slope that is 
greater than the corresponding one for real-space quantities {P~^), 
the more the bigger /3. This is expected since redshift distortions 
enhance linear density contrasts (eq.|17|l, and the amplitude of dis- 
tortions scales linearly with /3. We can therefore consider the factor 
/3~^ as the one coming from real-space dynamics, and the addi- 
tional factors involving /3 as redshift-space corrections. 

As stated above, in redshift space, unlike in real space, the 
velocity-gravity relation is not deterministic even in the linear the- 
ory, so aj^ > au > ap- However, the linear gravity and veloc- 
ity fields are tightly correlated, so Equations i24t-l|26|l yield val- 
ues of the ratio between and Vs very close to each other: e.g., 
for 13 = 0.512 (the value used in our numerical experiments, §|4) 
07^ = 2.579, au = 2.565, and ap = 2.552. 

To illustrate this point furthermore, one can examine the scat- 
ter on the conditional average {g)\v (mean g given v), where we 
have dropped out the subscripts 's' for simplicity of notation. For 
g and v being Gaussian-distributed, {g)\v ~ {vg)v / {v^), in agree- 
ment with linear regression, equation l[7]l. The scatter in this re- 
lation, (Tg|^, is the square root of the conditional variance, {(g — 
{g)\vY')\v The conditional variance is then equal to (1 — 
(see, for example, Appendix B of Ciecieljg et al. 2003). Thus, for 
Gaussian velocity and gravity fields in redshift space we have 



(1 - r>l 



(27) 



where r is the cross-correlation coefficient defined in eq. dllt . Note 
that the scatter is just one number, i.e., it is independent of the value 
of V. Equations (I21t-(|23|l yield 
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(28) 



In particular, for /3 = 1, r ~ 0.987, and for /3 = 0.5, r ~ 0.995: 
again, redshift-space velocity and gravity in the linear regime are 
(though not simply mutually proportional, like in real space) very 
tightly correlated. Inserting equation i28t into l l27t yields 



2 



/32 



12 
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(29) 



Hence, for (3 = 1, Gg\^ ~ 0.16CTg, and for /3 = 0.5, ~ 
O.lOag. This implies that the 'signal to noise ratio', S/N, of the 
estimate of a single galaxy's peculiar velocity from its gravitational 
acceleration can be as high as 10 {S/N ~ g/o'g\v ~ '^a/'^g\v — 
10). This is to be contrasted with the signal to noise ratio of the 
estimate of a galaxy's peculiar velocity from its distance and red- 
shift, which is typically below unity: an expected typical 20% rel- 
ative error in distanccQfor a galaxy at a distance of 40 Mpc0 
translates to the velocity error of 800 km ■ s~^, greater than typical 
peculiar velocities of galaxies. As a corollary, the intrinsic scatter 
in the redshift-space linear velocity-gravity relation is negligible 



^ See, for instance, Strauss & Willick (1995) for a review on issues related 
to peculiar velocity estimates. 

^ Where h = Hq/WO and Hq is the Hubble constant expressed in 
km/s/Mpc. 



compared to that introduced by the errors of measurements of pe- 
culiar velocities. 



3 VALIDITY OF THE DISTANT OBSERVER LIMIT IN 
THE REST FRAME OF CMB 

While equation ^ was written in the CMB rest frame, the choice 
of the reference frame for computing redshifts of galaxies can be 
more general. 



s = r + \v(r) 



(30) 



where Uorig is the (angle dependent) radial velocity of the origin of 
the system of coordinates (in the CMB rest frame). 
The Jacobian of eq. i30l yields (Kaiser 1987) 



1+5, 



(a)/ 



1 + 



v{r) — Vo 



1 + 



dv 
dr 



[l + ^'^^j .(31) 



This equation is valid as long as the mapping BOb does not induce 
any shell crossing. In general, non trivial singularities can appear in 
(5s*' (s), even if 5'®-* is finite. However, in practice, since one always 
performs some additional smoothing to the data, Si^\s) remains 
finite. The small r limit can be still problematic as it corresponds to 
a singularity in the system of coordinates. If the field v{v) is smooth 
enough, then in the neighborhood of the origin, one can write 



1 



v(v) — -Uorig 



1 + 



V(0) - Vo 



(32) 



Choosing a spherical coordinate system such that the 2 axis is par- 
allel to v(0) — Vorig, the term {1 + [v{r) — Uorig]/?"}"^ in eq. | |3U 
will create a singular surface of equation 



-|v(0) - Vorig|cos(6'), 61^72 



(33) 



that will concentrate to the origin, s = 0, in redshift space. Simi- 
larly the term (1 + dv/dr)^^ might become singular, but we see 
here that the situation nearby the origin is not different from what 
happens far away from it. 

The singular behavior of the form l l33t is expected to occur 
only in the neighborhood of the observer, but might be problem- 
atic when estimating the gravitational acceleration. While insert- 
ing this singularity in eq. lO, we notice that it should be an issue 
only for s = 0, where it coincides with the Green function sin- 
gularity (this is due to the finite mass of the singular surface). In 
other words, even though its small s behavior is difficult to pre- 
dict, the redshift space gravity should not be significantly affected 
by such a singular behavior at distances sufficiently large from the 
observer. In the CMB rest frame, Uorig ~ and v{0) would cor- 
respond to the Local Group velocity, say, |v(0)| ~ 630 km/s (see 
e.g., Erdogdu et al. 2006), so we would need r and therefore s large 
enough compared to 6 Mpc. More specifically, let us consider 
the real space sphere of radius r = |v(0)|. Its content is embedded 
in a sphere of radius s — 2|v(0) | in redshift space. The mass inside 
such a volume remains finite, but its internal distribution affects the 
gravity field at larger s in a non trivial way. What matters is that the 
multipole contributions of higher order than the monopole (the sub- 
structures within this volume) have negligible contribution on the 
gravity field. Using the wisdom from treecode simulations tech- 
niques (e.g. Barnes & Hut 1986; Barnes & Hut 1989), this amounts 
to 



2|v(0)| 



< 1. 



(34) 
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According to that criterion, the effect of the singular behavior near 
the origin should be of little consequence if s > 2|v(0)| = 13h^^ 
Mpc. 

If it is supposed now that s is indeed large enough, one can 
linearly expand eq. l |31b to obtain (Kaiser 1987) 



r or 



(35) 



(Note that, at linear order, one can assume s ~ r.) The distant ob- 



server limit would consist in dropping the term 2[v{r) 



from this equation. However, as extensively discussed in Kaiser 
(1987) and Kaiser & Lahav (1998), this term is in fact non neg- 
ligible in the redshift space gravitational acceleration as it induces 
the so-called rocket effect, resulting in a large r logarithmic diver- 
gence if Worig 7^ 0. This justifies the choice of CMB rest frame 
coordinate system, Uorig = 0. Still the remaining contribution of 
v{r)/r, although zero in average, might introduce some significant 
fluctuations on the large scale redshift space gravity field. In the 
linear perturbation theory framework, «(r) does not correlate with 
either (5*^^^ or dv/dr. As a result, while computing the sum of the 
fluctuations of 5^®' in a sphere of radius r = s' in eq. (|5j, what mat- 
ters is to see whether the fluctuations added by 2v{r)/r are small, 
or not, compared to the fluctuations added by 5^^^ — dv/dr. 
So let us estimate the ratio 



7^■ 



{[2v{v)/Tf) 



(36) 



{[5<^a) - dv/drf)' 

Calculations are similar to § [21 and one simply finds, in the linear 
regime. 



3 (1-f 2/3/3 -f /3V5) r2 P{k)k^dk ' 



(37) 



One can furthermore assume that the linear fields are smoothed, e.g. 
with a Gaussian window of size £. In that case one has to replace 
P{k) with P{k) exp[~{k£)^] in eq. ll37tFI For scale-free power 
spectra, P{k) oc fc", one finds using eq. (4.10b) of Bardeen et al. 
(1986) 



3(n + 3)(l + 2/3/3 + /3V5) 
which gives 

0.48 



7^■ 



0.48 

n + 3 \r) 



for p = 0.5. 



(38) 



(39) 



This means that smoothing increases the relative contribution of the 
term u(r)/r! In other words, smoothing makes the distant observer 
approximation worse. 

The standard cold dark matter (CDM) cosmology considered 
in the forthcoming numerical analyses assumes the non-relativistic 
matter density parameter Q, — 0.3, the cosmological constant A — 
0.7 and the Hubble constant i/o = 70 = 100 hkm-s'^ ■ Mpc"^ 
Using for instance the package of Eisenstein & Hu (1998) to com- 
pute P{k) with these cosmological parameters (and fJbaryons = 
0.04), one obtains numerically, in the absence of bias, 6=1, 



^ 0.17/i~^Mpc 

7c ~ , no smoothing, 



7^ 



2.7 h-^ Mpc 



lh~^ Mpc 



(40) 



(41) 



^ We suppose here that smoothing on v{r)/r is approximately equivalent 
to smoothing on f (r), prior to dividing by r, which should be reasonable. 



7^ 



10/i~^Mpc 



I = 10 /i"^ Mpc. 



(42) 



Clearly, for a catalog depth of the order of 70 Mpc as consid- 
ered below, we expect significant deviations from our theoretical 
predictions if the smoothing scale is as large as 10 Mpc. 



4 NUMERICAL EXPERIMENTS 

In this section we perform controlled numerical experiments to test 
the velocity-gravity relation, both in real and redshift space, on 
the dark matter distribution. These analyzes extend the work of 
Cieciel^g et al. (2003) who performed similar work but only in 
real space and on simulations using a pure hydrodynamic code ap- 
proximating the dynamics of dark matter, namely the Cosmolog- 
ical Pressureless Parabolic Advection code of Kudlicki, Plewa & 
Rozyczka (1996, see also Kudlicki et al. 2000). 

This section is organized as follows. In § 14.11 we present a 
new algorithm to interpolate the velocity and the gravity fields on 
a grid using an adaptive interpolating procedure inspired closely 
from smooth particle hydrodynamics (SPH). In i; 14.21 we describe 
the A''-body simulation set used in this work. In § 14.31 we discuss 
measurements on our dark matter samples and diluted counterparts, 
while mock galaxy catalogs are addressed in §|4!4] 

Since we assume flat cosmology, + A = 1, we adopt the 
approximation 



(43) 



for function fip), as it is known to be slightly more accurate in 
that case than the traditional fit (Bouchet et al. 1995). Finally, 
let us recall that 'DOL' means "distant observer limit". 



4.1 Algorithm used to estimate velocity and gravity fields 

In real observations, the catalogs used to estimate cosmic velocity 
and gravity are in general different. In particular, the velocity field 
tracers are not necessarily representative of the underlying density 
field. Here we assume to simplify that the same catalog is used 
to probe velocity and gravity fields and that all the objects in the 
catalog have the same weight, or from the dynamical point of view, 
the same mass. Also, since gravitational force is of long range, it is 
necessary to estimate it in a domain Vgrav large enough compared 
to the effective volume V where the velocity-gravity comparison 
is actually performed. We assume here that it is indeed the case. In 
this paper, Vgrav will simply be the simulation cube, while V will 
be a sub-volume included in this cube. 

The basics of our method is to perform adaptive smoothing on 
the particle distribution to obtain a smooth velocity and a smooth 
density field on a regular grid encompassing Vgrav and of resolu- 
tion A'g, while preserving as much as possible all the information. 
Then additional smoothing with a fixed window, preferably Gaus- 
sian, can be performed a posteriori in order, e.g., to be in the lin- 
ear regime, since our analytic predictions in principle apply to that 
regime. We shall see however below that this additional step is not 
necessarily needed and can in fact complicate the analyzes (see also 
§[3}. 

Whether we perform the measurements in real or redshift 
space does not change our approach: we assume that we have at 
our disposal a set of TV points with three-dimensional coordinates 
and a scalar velocity. This latter is either the radial velocity when 
we consider redshift space measurements, or the z coordinate of 
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the velocity vector when we consider real space measurements or 
redshift space measurements in the DOL approximation. 

The main difficulty is to reconstruct a smooth velocity field 
on the sampling grid. Let us remind as well that the velocity field 
we aim to estimate is a purely Eulerian quantity, in other words, a 
mean flow velocity. In particular, what we aim to measure, in terms 
of dynamics, is something as close as possible to a moment of the 
density distribution function /(fl, vl) in phase-space: 



v(r) = ^ / d''vLVL/(r,VL) 



(44) 



where p is density, the source term of the Poisson equation to esti- 
mate the gravitational potential, given by 



P(r) 



d VL/(r, vl) 



(45) 



Note that equation ( I44l > applies only to real space, but we shall 
come back to that below. 

Since we aim to estimate velocity and density on a finite res- 
olution grid, a more sensible way of performing the calculations 
is to integrate equations i44t and l l45t on a small cubic patch Ar'' 
corresponding to a grid element to reduce noise as much as pos- 
sible and to make the calculation conservative (i.e., the total mass 
and momentum is conserved). In practice, we do not perform this 
calculation exactly, but only approximately using an approach in- 
spired from smooth particle hydrodynamics (SPH, see, e.g. Mon- 
aghan 1992), as we now describe in details. 

In the SPH approach, each particle is represented as a smooth 
cloud of finite varying size 2i?sPH depending on local density, i.e. 
on the typical distance between the particle and its A'^sph closest 
neighbors (which can be found quickly with e.g. standard KD-tree 
algorithm), where A'^sph is usually of the order of a few tens. As a 
result, smooth representation of the density and velocity fields can 
be obtained at any point of space by summing up locally the contri- 
bution of all the clouds associated to each particle. With appropriate 
choice of the SPH kernel, these functions can be easily integrated 
over each cell of the grid. 

The problem with this approach is that it does not guaran- 
tee that local reconstructed density is strictly positive everywhere, 
which can leave regions where the velocity field is undefined. To 
solve that problem, we start from the grid points, which we assume 
to be virtual particles for which we find the A'^sph closest neighbors 
to define the SPH kernel associated to this grid site. However, one 
has to take into account the fact that the sought estimates should 
roughly correspond to an integral in a small cubic patch Ar"^. In 
particular, all particles belonging to a grid site should participate 
to such an integral. This issue can be addressed in an approximate 
way as follows: 

(i) Count and store for each grid site, the number A^p of particles 
it contains; 

(ii) If A^p ^ A'^sPH, then perform SPH interpolation at the grid 
site as explained below using the A^'p particles contributing to it 
instead of the A^'sph closest neighbors; 

(iii) If < A'^sPH, then find for the grid site the A'^sph closest 
particles and perform SPH interpolation at the grid site as explained 
below. 

Given a choice of the SPH kernel, S{x) [which should be a 
monotonic function verifying <S(0) — 1 and S{x ^ 2) = 0], and 
a number A'x of neighbors, with TVx = A^p or A'^sph according to 
the procedure described above, the interpolation of a quantity A on 
the grid site is given by 



A{i,j,k) = 



JVx-i 

E 



AlWiS 



d. 



Rspii{i,j,k) 



.(46) 



In this equation, Ai is the value of A associated to each particle 
^SPH (* , j, k) is half the distance to the furthest neighbor of the grid 
site k) among the TVx, d; is the distance of the Z-th particle to 
the grid site, and Wi is a weight given to each particle such that 
the total contribution of every particle to all the grid sites is exactly 
unity: 

1 di 



Wi = 1/Si 



E 



[i?SPH(i, j, \ 7?sPH(i,i, fc) 



In practice, the interpolated density thus reads 

p{i,j,k) = m(i, j,fc)/Ar^, 



(47) 



(48) 



by taking Ai — mi in equation l l46b . where mi is the mass of each 
particle. An interpolated velocity reads 



v{i,j,k) =p{i,j, k)/m{i,j, k), 



(49) 



where p{i, j, k) is the interpolated momentum in cell (i, j, k). It is 
obtained by taking Ai = pi = miVi in equation {46), where vi is 
the velocity (the radial component, or a coordinate) of particle I. 

Note that there can be some particles for which Si = in 
equation ( 147 1 (in that case Wi = by definition), i.e. which do 
not contribute at all to the interpolation. In the practical measure- 
ments described later, this can happen at most only for a very small 
fraction of the particles, typically of the order of 0.1 percent, and 
affects the results insignificantly. There is, however, another notice- 
able defect in our method, due to the unsmooth transition between 
the two interpolation schemes when the number of particles A^p per 
grid site becomes larger than TVsph, that affects in a non trivial way 
the interpolated density, but again, it does not have any significant 
consequences for the present work. 

In practice, we take A'sph = 32 for all the measurements 
described below and the following spline for the SPH kernel (Mon- 
aghan 1992) 

S{x) = 1 - 1.5 X- + 0.75 x^ OsCxsJl 
S{x) = 0.25 (2 - x)^ 1 s; a; ^ 2 

S{x) = 0, x^2. (50) 

Once the density field is interpolated on the grid, it is easy 
to estimate the gravitational potential from it by solving the Pois- 
son equation in Fourier space, keeping in mind that the edge of the 
volume V, where the velocity-gravity relation is tested, should be 
sufficiently far away from the edges of the volume Vgiav, which 
is itself included in the grid. It would be beyond the scope of this 
paper to discuss other problems related to incompleteness or edge 
effects. The main one is related to the uncertainties on the gravita- 
tional potential induced by the obscuring due to our own Galaxy. 
While this can be certainly an issue, other incompleteness prob- 
lems such as segregation in luminosity can be addressed by giving 
the proper weight (or mass) to the galaxies in the catalog. This of 
course needs strong assumption on the bias, and can work only for 
populations of not too bright galaxies. 

Note finally that, as explained before, we perform "naively" 
our SPH interpolation whether we are working in real space or 
in redshift space. The fact that we use SPH interpolation in red- 
shift space is sufficient for our purpose as long as we are in the 
DOL limit, although it does not correspond anymore exactly to a 
simple moment of the Vlasov equation. However, the interpolation 
becomes somewhat questionable when the assumption of DOL is 
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dropped: the nature of the interpolation changes with distance from 
the observer. In particular, if we assume that the smoothing ker- 
nel was fixed, projections, e.g. passing from the velocity vector to a 
one dimensional quantity such as the radial velocity, and smoothing 
do not commute anymore. The same problem arises for commuta- 
tion between the calculation of the radial part of the gravitational 
force and smoothing. It is therefore necessary to carefully check 
that the simplistic nature of our interpolation does not introduce 
any systematic bias in the measurements, when performing them in 
redshift space relaxing the approximation of a distant observer. 



4.2 The A'^-body simulation set 

We performed a high resolution simulation using the adaptive 
mesh refinement (AMR) code RAMSES (Teyssier 2002). As al- 
ready mentioned earlier, the cosmology considered here assumes 
n = 0.3, A = 0.7andi?o = 70 = 100 /i km- ■Mpc~\ Initial 
conditions were set up using Zerdovich approximation (Zel'dovich 
1970) to perturb a set of particles disposed on a regular grid pattern 
to generate initial Gaussian fluctuations with a standard ACDM 
power spectrum. To do that, we used the COSMICS package of 
Berstchinger (1995). The simulation involves 512"^ dark matter par- 
ticles on the AMR grid initially regular of size 512'^, in a periodic 
cube of size 200 Mpc. Then, additional refinement is allowed 
during runtime: if cells contain more than 40 particles they are di- 
vided using standard AMR technique (with a maximum of 7 levels 
of refinement). Note finally for completeness that the normaliza- 
tion of the amplitude of initial fluctuations was chosen such that the 
variance of the density fluctuations in a sphere of radius 8h^^ Mpc 
extrapolated linearly to the present time was given by as = 0.93. 



4.3 Measurements on the dark matter distribution 

From our RAMSES simulation, we extracted a number of dark- 
matter samples, as described in details in Table[T] namely: 

(i) A high resolution grid sample of size A'g = 512 for testing 
the velocity-gravity relation in real space; 

(ii) The analog of (i) for testing the velocity-gravity relation in 
redshift space, using the DOL approximation; 

(iii) The analog of (ii), but without using the DOL, for testing the 
validity of this approximation. In that case the radial coordinate of 
the velocity and gravity field was estimated using an observer at 
rest at the center of the simulation box. The comparison of veloc- 
ity to gravity was performed in a sphere of radius 70 h^^ Mpc to 
avoid edge effects due to the loss of periodicity while projecting in 
redshift space; 

(iv) A set of 125 low resolution grids (A'^g = 128) in redshift 
space, for estimating (at least partly) cosmic variance effects and 
effects of structures nearby the observer. These samples were gen- 
erated by locating the observer in the simulation box on a regular 
grid of size 5x5x5. Again, in these samples, the velocity-gravity 
relation was tested in a volume of radius 70 Mpc centered on 
the observer (exploiting the periodic nature of the simulation box). 

Since the catalogs we consider in point (iv) represent a significant 
fraction of the simulation volume, we know however that effects of 
cosmic variance are likely to be underestimated. 

We now perform a visual inspection of the fields, followed by 
measurements of the joint PDF of velocity and gravity, first in real 
space, then in redshift space. 



4.3.1 Visual inspection 

Figure[T]shows the density, the gravity field and the velocity field in 
a thin slice extracted from the simulation, both in real space and in 
the DOL redshift space. Despite our interpolation procedure, there 
are some minor discreteness artifacts left, visible on the density 
field. One can notice for example a few underdense regions where 
the initial grid pattern, distorted by large scale dynamics, is still 
present. Such artifacts do not show up on the smoother z coordinate 
of the gravity field. On the other hand, this latter seems to suffer 
from a few aliasing defects in real space, as vertical spurious lines 
in the vicinity of deep potential wells. We did not bother trying 
to understand such aliasing effects, because they have negligible 
impact on the measurements^ Interesting to notice here is the nice 
agreement between the z coordinate of the velocity field and the z 
coordinate of the gravity field in underdense regions. These regions 
dominate the velocity-gravity statistics, since we have a volume 
weighted approach. Note the particular features in the velocity field 
associated to the filaments in the density distribution, as well as the 
'fingers of God' (FOG) in redshift space. These FOG, which are 
mainly associated to the dark matter halos, are expected to induce 
some particular properties on redshift space statistics: 

• First, and this is a straightforward consequence of passing 
from real space to redshift space, there is, on the velocity field, 
an 'inversion' effect inside FOG. In other words, inside a finger of 
God, the variations of the velocity field are opposite to what hap- 
pens in its nearby environment. This can be explained in the follow- 
ing way. Inside a halo, which can be, in real space, considered as a 
point-like structure in first approximation, there are particles with 
positive velocities and particles with negative velocities. Assume, 
to simplify, that this halo has itself a zero center of mass velocity. 
It will, in redshift space, look like an elongated structure. Particles 
in this structure that have positive z velocity will be above the cen- 
ter of mass, and particles with negative velocity will be below the 
center of mass. As a result, one expects, for the z coordinate of the 
velocity, the finger of God corresponding to this halo to have posi- 
tive velocities (towards light color on bottom right panel of Fig.lTJ 
above the center of the halo and to have negative velocities (towards 
dark color) below the center of the halo. This is indeed what we can 
observe. However, in the nearby environment of the halo, the situ- 
ation is somewhat opposite. Indeed, the halo is expected to lie in a 
filament, which itself represents a local potential well. If this well 
is not too strong to induce shell-crossing while passing from real to 
redshift space, the general trend of the velocity field is not modified, 
compared to real space: negative sign (towards dark color) above 
the filament, positive sign below it. If the filament corresponds to a 
sufficiently deep potential well such that shell crossing occurs, then 
the same effect than for FOG is expected, as can be noticed on the 
bottom right panel of Fig.[T]for the largest filaments. 

• Second, due to the FOG stretching, halos are more elongated 
in redshift space, and a natural consequence is that the correspond- 
ing potential well is less deep. That is why the z coordinate of the 

* These aliasing effects might be related to some minor defect in the pai- 
ticular Fourier transform algorithm we are using (Teuler 1999), but this hy- 
pothesis seems to be contradicted by some accuracy tests performed on this 
algorithm (Chergui 2000). A more sensible explanation is that these effects 
ai'e induced by the way we interpolate the density field, in particular by the 
transition between the A^p ^ A^spjj and the A'^p < Afgpjj regimes, com- 
bined with the fact that our Fourier Green function for computing the grav- 
itational acceleration is simply proportional to k/fc^, without additional 
filtering. 
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Figure 1. The density (upper panels), gravity (middle panels) and velocity field (lower panels) in a one-mesh-element-thick slice extracted from the samples 
(i) and (ii), which correspond to a grid resolution of 512^. The left and right columns of panels correspond respectively to real space and distant observer 
redshift space. In the top panels, the color scale is logarithmic, from white to black. In the foiu" bottom panels, the z coordinate of each field is displayed. The 
color scale is linear and normalized in such a way that direct visual comparison between the four panels can be performed. One can notice clearly 'fingers of 
God' on the upper right and lower right panels. In such 'fingers of God', the velocity field tends to have a different signature compared to that in the nearby 
environment, as explained in the text. Finger-of-God effects also tend to soften the gravitational potential at small scales (cf. middle panels). 
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Table 1. Summary of the characteristics of the dark matter samples used in this paper and the con'esponding main quantitative results. The first column gives 
the name of the sample as used in the text. The second one mentions whether the measurements are performed in real or redshift space, the acronym DOL 
meaning that the approximation of an infinitely remote observer was used. The third column gives the number A^rea of realizations considered. The fourth 
column gives the size of the sample, which is the full periodic simulation cube of size L when the letter L is used and a sphere of radius R when the letter 
R is used. Scales are expressed in units of Mpc. The fifth column gives the resolution of the rectangular grid (matching the simulation box) used to 
sample various fields from the particle distribution. The sixth column gives the number A^obj of objects that are actually contained in the sample. A star means 
that this number is only approximate, con'esponding to a simple rescaling of the total number of objects in the periodic cube by taking into account the size of 
the sampled volume, a sphere of radius R.. The 7th, 8th and 9th columns give the measured SI using {gf) / (viQi) (1st number), {9i) / (^^) (2nd number) 
and (viQi) I (v^ ) (3rd number) as estimators of the velocity-gravity relation, where i = s if the measurements are performed in redshift space and i = zii the 
measurements are performed in real space or distant observer redshift space. In principle, all these estimators should give the same answer (nearly the same 
answer in redshift space) in the ideal case where the linear theory applies, and various sources of noise (such as shot noise, cosmic variance, systematics, etc.) 
are negligible. The disagreement between these various estimates can be used as a proxy to estimate error bars on actual measurements. Where possible, an 
error bar obtained from the dispersion among various reahzations is quoted as well. The last column explains what exactly the quoted numbers correspond to, 
in particular if the fields were smoothed with a Gaussian window of size £ = 10 h^^ Mpc ("smoothed") or not ("no smoothing"), if all the PDF was used to 
perform the measurements ("all") or only the most likely region defined by Equation i52\ ("1.5a" isocontour"). The symbol (■ ■ ■) means that an average was 
performed over the number of available realizations, when relevant. 



Samp. 


Content A^rca 


Size 






OOJ 


Measured Q from: 

(vg) y {v'^) (ii^> 


Comment 


(i) 


real space 1 


L = 


200 


512 


5123 


0.103 


0.157 


0.240 


all, no smoothing 














0.251 


0.257 


0.264 


all, smoothed 














0.299 


0.303 


0.307 


1.5(T isocontour, no smoothing 














0.295 


0.297 


0.299 


1.5(T isocontour, smoothed 


(ii) 


DOL redshift space 1 


L = 


200 


512 


5123 


0.143 


0.280 


0.578 


all, no smoothing 














0.206 


0.223 


0.242 


all, smoothed 














0.296 


0.303 


0.311 


1.5(7 isocontour, no smoothing 














0.274 


0.279 


0.284 


1.5(7 isocontour, smoothed 


(iii) 


redshift space 1 


R = 


70 


512 


2.41e7* 


0.098 


0.287 


0.987 


all, no smoothing 














0.167 


0.197 


0.231 


all, smoothed 














0.316 


0.332 


0.349 


1.5(7 isocontour, no smoothing 














0.269 


0.284 


0.300 


1.5(7 isocontour, smoothed 


(iv) 


redshift space 125 


R = 


70 


128 


2.41e7* 


0.115 


0.280 


0.817 


( all ) , no smoothing 














0.023 


0.030 


0.269 


en'or from dispersion 














0.177 


0.200 


0.226 


( all ) , smoothed 














0.022 


0.021 


0.023 


error from dispersion 














0.283 


0.294 


0.304 


( 1.5(7 isocontours ), no smoothing 














0.039 


0.040 


0.041 


en'or from dispersion 














0.255 


0.262 


0.270 


( 1.5(7 isocontours ), smoothed 














0.031 


0.032 


0.034 


error from dispersion 



gravitational force is less contrasted in the middle right panel of 
Fig-dthan in the middle left one. 

Still, we can see that both in real and in redshift space, the velocity- 
gravity relation is going to be in quite good agreement with ex- 
pectation from linear theory. Indeed and again, our measurements 
are volume weighted, dominated by underdense regions. The color 
scale on four bottom panels of Fig.[T]has been chosen such that if 
linear theory applies, the same color pattern should be found for the 
gravity and velocity field, which is the case at first glimpse, except 
in the densest regions, corresponding to halos or rich filaments. 

Finally, similarly as the right column of panels of Fig.[T] Fig.[2] 
displays the density, the gravity field and the velocity field in red- 
shift space, but without using the DOL. If one now takes into ac- 
count the radial nature of the projection, it is clear that the conclu- 
sions of the previous discussion remain unchanged, at least at the 
qualitative level. 



4.3.2 Quantitative comparison in real space 

Figure |3] shows in grey scale the measured joint probability distri- 
bution function (PDF) of the z coordinate of gravity and velocity 
fields extracted from the real space sample (i). The striking result 



is that the regions of the best likelihood (larger values of the PDF, 
darker places) match very well the prediction given by linear the- 
ory (thick solid line), even in the highly non-linear regime (upper 
panel). This remarkable property is mainly related to our volume- 
weighted approach: results are mainly influenced by underdense 
regions, which are weakly evolved from the dynamical point of 
view and are expected to match well linear theory predictions. This, 
and the 'propeller' shape of the bivariate PDF, can be understood 
in more detail by using the spherical top-hat model as a proxy of 
nonlinear dynamics. In that case, up to shell crossing, the velocity- 
gravity relation reads approximately 



R 
3 



PR 



3/2' 



(51) 



(this can be easily derived from Bemardeau 1992; 1994), where R 
is the distance from the center of a top-hat fluctuation (in km ■ s^^). 
This equation is valid inside the fluctuation, which can be over- 
dense (negative v) or underdense (positive v). One can then imag- 
ine, to simplify, the density field as a patchwork of spherical fluctu- 
ations, which correspond to a set of curves given by Equation ( 151b . 
as shown in the top panel of Fig. [3] Note here that we should, 
for the picture to be correct, take into account the fact that we 
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Velocity 

RAMSES, Real space, Lg=0.4h"' Mpo 




Velocity 

Figure 3. The velocity-gravity relation in real space as a scatter plot, compared to theoretical prediction in the linear regime (thick diagonal line), as measured 
from sample (i) extracted from the simulation using the z coordinate of each field. The color scale is defined as follows: the darkest area is encompassed by a 
68 percent isocontour (eq. |52) . corresponding approximately to a l.Scr contour in the Gaussian case. Then, lower isocontours, corresponding to lighter color, 
scale logarithmically, which emphasizes the tails of the PDF. 

The two upper panels con'espond to the raw data interpolated on the grid with the SPH-like interpolation detailed in § |4. 1| while in the bottom panel, an 
additional Gaussian smoothing with kernel of a fixed size t = lOh^^ Mpc was peii'ormed. The smoothing scale is indicated on three panels (in the two top 
panels, it corresponds to the mesh element size). In the top panel, predictions of velocity-gravity relation according to the spherical top hat model are displayed 
as solid curves. As discussed in detail in the text, these curves help to understand why the linear prediction still represents the most prominent feature of the 
relation, even in the highly nonlinear regime, while the tails due to nonlinear dynamics induce a 'propeller' shape of the scatter plot, as already noticed by 
Ciecieljg et al. (2003). Such a propeller shape results in a bias on the measurement of the slope by using directly second-order moments of the bivariate PDF. 
This is demonstrated by plotting on middle (and bottom) panel the dotted line which gives the estimated slope from the ratio -y/ (ffl)/ (^D- This dotted line 
should be compared to the thick solid line. The dashed lines correspond to slopes obtained by the conditional averages (9z)/((?zf z) and (g^ j;z)/(d^). In the 
middle panel these slopes are significantly different, consistent with significant amount of the scatter in the relation, due to nonlinear effects. In the middle 
panel, also the 1.5(T contour expected in the Gaussian case is shown. 
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are using only the z coordinate of the fields, Qz = g cos 6 and 
= V cos 8, where 6 is the angle between the z-axis and the ra- 
dial vector R. In the top panel of Fig. |3] we consider the cases 
7?cos6l = ±10, ±100, ±1000, ±3000 km/s. Each curve has a 
stopping point corresponding to the maximum possible value of v, 
■Umax = I3R/2, which reflects the fact that there is an expected up- 
per bound for the expansion speed of voids. This property excludes 
the upper left and lower right quadrants of the velocity-gravity di- 
agram to be populated too far from linear theory prediction. On 
the other hand, overdense fluctuations tend to populate the upper 
right and the lower left quadrants, above and below linear the- 
ory prediction, respectively. Furthermore, since the low-u regime 
converges to linear theory, all the curves corresponding to Equa- 
tion ( 1511 ) superpose in that regime, creating a 'caustic' of best like- 
lihood nearby the maximum of the joint PDF, explaining the very 
good agreement with linear expectation in that region. As a result, 
we now understand, thanks to the spherical top-hat model, both the 
'propeller' shape of the bivariate PDF, as well as the remarkable 
agreement with linear theory prediction nearby its maximum, even 
in the highly non-linear regime (see also Cieciel^g et al., 2003). 
The arguments developed here are oversimplified, but capture the 
main features of the dynamics of the large scale structures prior to 
shell crossing in real space. Beyond shell crossing, there is a mix- 
ing effect that tends to decorrelate velocity from gravity, implying 
a widening of the bivariate PDF and even larger tails in the upper 
right and the lower left quadrants. However such an effect does not 
affect significantly the region of best likelihood. 

A straightforward consequence of the above discussion is that 
the joint PDF of gravity and velocity is substantially non Gaussian 
due to non-linear contributions in the dynamics, as supported by 
the examination on middle panel of Fig. [3] of the 1. 5a elliptic iso- 
contour of the Gaussian distribution with the same second order 
moments as the measured PDF. As a result, the direct measurement 
of (3 parameter from the moments of the joint PDF using linear the- 
ory predictions is biased to lower values, due the propeller shape 
of the PDF. This is illustrated on middle panel of Fig.[3]by the dot- 
ted line, which gives the velocity-gravity relation obtained from 
the second order moments of the PDF, {g'i)/{'vl), while the two 
dashed lines correspond to the conditional averages {g1) / {gzVz) 
and {gzVz) / {v1). In the linear regime, let us remind from §|2]that 
these three quantities should be equal. Due to non-linear effects, the 
relation between gravity and velocity is not deterministic anymore, 
kaDlkg^Vz) > {gi)/(vi) > (ff^w^)/(uf>, and /? is underesti- 
mated, leading to an effective bias or an underestimated value of f2, 
as shown in Table [T] Note, as expected, that additional smoothing 
of the fields with a Gaussian window of radius 10 h^^ Mpc helps 
to reduce non Gaussian features as well as this effective bias, and 
makes the overall relation between gravity and velocity tighter. 

However, Fig.[3]appeals for a more sophisticated way to mea- 
sure the velocity-gravity relation than simply using directly the 
moments of the joint PDF. Since the region around the maximum 
seems to agree rather well with the linear prediction, this suggests 
to estimate the moments from the PDF only within that region. To 
perform such an exercise, we selected the 68-percent PDF isocon- 
tour, Ves, such that 

/ dgdvP{g,v)= 0.68. (52) 

Jv{g,v)>Ve8 

This corresponds roughly to a l.Scr contour in the Gaussian case. 
The PDF is then set to zero outside the contour of best likeli- 
hood. From this truncated PDF, the moments are estimated again, 
leading to a much better estimate of fl, agreeing at the few per- 



cent level with the true value, as shown in Table [T] regardless 
whether additional Gaussian smoothing with a 10 Mpc size 
window is performed on the fields, or not. Furthermore, the estima- 
tors ^yjgW/Jvf}, {gz)/(gzVz) and {gzVz)/{vz) differ now only 
slightly from each other, reflecting the narrowness of the region of 
best likelihood around the linear expectation. 

4.3.3 Quantitative comparison in redshift space 

Figure |4] is similar to Fig. [3] but shows the results obtained from 
the sample (ii), in the DOL redshift space. The same striking result 
obtained as in real space holds: the region of maximum likelihood 
agrees very well with the linear theory prediction, which gives a 
velocity-gravity slope larger than in real space because of the en- 
hancement of large scale density contrasts due to projection in red- 
shift space (symbol LSRD - Large Scale Redshift Distortion - on 
the left panel). 

However, there is on Fig. |4] a noticeable new feature on the 
PDF visible on left panel, in addition to the propeller shape: the 
joint PDF seems now to present tails in the directions orthogonal 
to the maximum likelihood domain. Note that these tails tend to 
disappear with smoothing which then makes the bivariate PDF look 
very much like the real space one. 

This new feature is due to the finger-of-God effects already 
discussed at length in ij |4.3.1l These effects have mainly two conse- 
quences: first, the gravitational potential is less contrasted at small 
scales due to the elongated nature of fingers of God, which tend to 
reduce the amplitude of the gravitational acceleration inside them 
and in their neighborhood. Second, there is an inversion effect 
which changes the sign of the velocity. The consequence is that 
some points that were, in real space, on the upper right quadrant 
are moved to the upper left quadrant as indicated by the arrow on 
the left panel of Fig.|4] These points create a tail on the upper left 
quadrant; similarly some points that lied in the lower left quadrant 
contribute to the tail on the lower right quadrant. 

In fact, in the same way as we did in real space, the results 
obtained in redshift space can again be interpreted qualitatively us- 
ing the top hat spherical collapse model, as detailed in AppendixlAl 
Note interestingly that in the linear regime, the top hat model gives 

a value significantly smaller than what is obtained from Eqs. ( 124b . 
( 125 1 and l l26t : equation l l53b gives gz/vz — 2.3 instead of 2.6. This 
disagreement with the statistical expectation from linear theory is 
explained in Appendix [A] However, despite the limitations of the 
spherical top hat model, the arguments developed previously in the 
real space case to explain why the linear regime dominates the most 
likely part of the joint PDF still hold. 

Thanks to finger-of-God effects, if no additional smoothing 
is applied to the interpolated fields, the measured joint PDF is now 
much more symmetric about the linear prediction than in real space. 
As a consequence, the slope obtained from the direct measure- 
ment of (g'i) /{v'i) (dotted line on left panel of Fig. |4j agrees 
now well with linear theory, while we still have {gl) / {gzVz) > 
\J (ffz) / ("z) > (ffz'^z)/(^2)- However, when an additional Gaus- 
sian smoothing is performed, the fingers of God tend to be diluted. 
This implies on Fig. |4] a clockwise rotation of the upper left tail 
to the upper right and of the lower right tail to the lower left. As 
a consequence, with our 10 Mpc scale Gaussian smoothing, 
one converges to similar behavior as obtained in real space, with 
the propeller effect significantly biasing the overall slope, implying 
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RAMSES, Eedshift space DOL, Lj=0.4h ' Mpc RAMSES, Redshift space DOL, L=10h"' Mpc 




Velocity Velocity 

Figure 4. The velocity-gravity relation in distant observer redshift space as measured from sample (ii) extracted from the simulation, using the z coordinate 
of the fields. The color scaling is the same as in Fig. [5] Again the match between the most likely part of the bivariate PDF (darkest region) and the prediction 
from linear theory (thick solid line) is close to perfect. The dotted line gives the slope obtained directly from the ratio ^ (gf ) / (»;| ) . The left panel which 
coiTesponds to minimum amount of smoothing can be compared to the two upper panels of Fig. [3] The differences with real space results are explained by 
the arrows: large scale redshift distortion (LSRD) increases the effective slope of the velocity-gravity relation, while the finger of God (FOG) effect moves 
points from the upper right down to the left, and points from the lower left up to the right, creating two additional tails of the PDF. As discussed in the text, 
this is due to the inversion of the signature of velocity happening in fingers of God as seen in the bottom right panel of Fig.[T] and to the effective smoothing 
of gravitational potential due to the extended nature of the FOG. When smoothing with a Gaussian window of radius ^ = 10 /i~^Mpc is performed, these 
two tails move back to the upper right and the lower left of the panel: one recovers the propeller shape obtained in real space. However, its effect is more 
pronounced, because there is still some remnant of the anti-diagonal effect, depending on the level of smoothing. 



again that Q, is biased low. This bias is more pronounced in red- 
shift space than in real space, because there is still some remnant 
of the anti-diagonal effect, depending on the level of smoothing. Of 
course, smoothing at larger scales would make the agreement with 
linear theory prediction better again. 

To reduce the bias, one can play again the exercise of measur- 
ing the slope of the velocity-gravity relation by selecting the region 
of best likelihood, as in Equation ( I52l l. The agreement with linear 
theory prediction is improved as expected, and the corresponding 
measured value of matches very well the real value as illustrated 
by Table [T] However, additional smoothing tends to mix the non- 
linear finger of God effects with linear features, contaminating the 
region of best likelihood. As a result, the measured value of Q. is 
slightly biased to lower values (0.28 instead of 0.3) and changing 
the likelihood contour selection does not improve significantly the 
results. This mixing effect brought by smoothing suggests that in 
fact it might be not so wise to perform additional smoothing of 
the interpolated fields. Even though smoothing at sufficiently large 
scales brings better overall agreement with linear theory, it makes 
the measurements much more sensitive to finite volume and edge 
effects and it is furthermore not needed here since linear theory 
regime always dominates the region of best likelihood. Besides, we 
showed in §|3that smoothing was making the DOL approximatioi{3 
worse. 

Figure |5] is exactly the same as Fig. [4] but the measurements 
are now performed in sample (iii), without assuming that the ob- 
server is (infinitely) remote. Qualitatively, all the conclusions de- 
rived from the analysis of Fig. |4] still hold. The main difference is 
that the measurements are more noisy due to the size of the sam- 
pled volume, now 5 times smaller. Therefore, the recovered value 



^ Valid, however, by definition, for the sample currently under considera- 
tion, (ii). 



of Q. from the measurements of the moments in the region of best 
likelihood of the joint PDF of the unsmoothed fields is good only at 
the ten percent level (Jl = 0.33 instead of 0.3, see Table[TJ. At this 
level of accuracy, we find that the linear prediction, which was de- 
rived in the approximation of DOL, agrees well with the measure- 
ments. To improve the quality of this comparison. Figure |6] gives 
the same scatter plots as in Fig.|5] but for sample (iv), i.e. after co- 
adding the contributions of 125 different observers. The value of 
derived from the best likelihood region of left panel of Fig.[6]agrees 
now at a few percent level with the true value of f2 (the estimated 
value is 0.294, see TablelTJ, showing that the distant observer limit 
is an approximation good enough for deriving the linear prediction. 
Note again the bias to lower values (0.262) on the measured value 
of Vl brought by additional Gaussian smoothing with a 10 ft^^ Mpc 
Gaussian window. This bias is more pronounced than in the DOL 
sample (ii), and this is certainly at least partly due to the fact that, 
as discussed in § [S] deviations from the DOL limit are not negli- 
gible anymore for such a smoothing scale, given the sample depth, 
and they add to the mixing between linear and nonlinear features 
discussed above. It could also come partly from the edge effect dis- 
cussed in caption of Fig. [6] 

The dispersion among the 125 different observers leads to a 
typical error on of the order of 0.04, suggesting that the errors 
related to the choice of the position of the observer - which probes 
the space of configurations for the s = singularity discussed in 
§[3]~ are of the order of 10 percent for our i? = 70 h"^ Mpc ra- 
dius catalog. These errors include cosmic variance effects, but these 
latter are probably underestimated because our spherical redshift 
samples represent a rather significant fraction of the total simula- 
tion volume. 

Nevertheless these measurements illustrate the relative robust- 
ness in that respect of our velocity-gravity estimator to derive SI 
from a large-scale galaxy survey. However, we reiterate that, as al- 
ready mentioned in beginning of § 14.11 the volume used to com- 
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Velocity Velocity 

Figure 5. Same as in Fig.|4] but for the sample (iii) which does not assume that the observer is (infinitely) distant. Although the general behavior is exactly the 
same as in Fig.|4] several explainable differences can be noticed. The measurement is more noisy and the apparent width of the bivariate distribution is smaller 
than on Fig.|4] This is mainly related to the fact that the sampled volume covers only 18 percent of the simulation box. Still, the most likely part of the PDF 
fits very well the slope predicted by linear theory in the infinite remote observer approximation, lending credence to our simplified approach. 



RAMSES, Redshift space, Lj=1.6h"' Mpc RAMSES, Redshift space, L=10h"' Mpc 




Velocity Velocity 



Figure 6. Same as in Fig.|5] but we now consider sample (iv), where 125 different observer positions are located on a regular grid in the simulation box and 
where the interpolated field is computed only on a 128^ grid. As expected, the resulting bivariate PDF is smoother than in Fig. [5] since the simulation volume 
is now fully sampled by all these observer positions. 

Note the significant asymmetry of the tails of the bivariate PDF around the major axis of the maximum likelihood region. This is an expected consequence 
of edge effects. Indeed, while performing redshift space projection to construct a spherical catalog of finite radius iimax, the galaxies near the edge of the 
catalog with positive radial velocity and radial position r < i?max tend to get outside the sample. On the other hand, galaxies near the edge of the catalog 
with negative radial velocity and radial position r > Rmax tend to get inside the catalog. As a result one expects more galaxies with negative radial velocities 
than with positive ones, hence the asymmetry of the PDF. This effect mainly affects fingers of God, so the anti-diagonal tails of the bivariate PDF. Importantly, 
this effect does not affect too much the region of best likelihood. 



pute gravity field should be significantly larger than the actual vol- 
ume used to perform the velocity-gravity comparison. The present 
analyses suppose it is the case as indeed achieved by most recent 
three dimensional galaxy catalogs such as the 2MRS, for which 
the median redshift corresponds to a half-depth of ~ 6Qh~^ Mpc 
(Erdogdu et al. 2006). 



4.3.4 Ejfects of dilution 

Up to now, we have measured the velocity-gravity relation in very 
rich catalogs, where the number of objects was so large that dis- 
creteness effects could be considered as negligible. In real galaxy 



catalogs, the number of objects is much smaller, particularly when 
tracers of the velocity field are taken at concern. Before addressing 
the issue of biasing between the galaxy distribution and the dark 
matter distribution, we therefore examine pure discreteness effects, 
by diluting our dark matter samples. This dilution will not only 
bring a shot noise contribution, it will also increase the overall level 
of smoothing which is performed by our interpolation procedure. 

To be able to quantify as accurately as possible the biases in- 
duced by discreteness in the mock "galaxy" catalogs considered in 
next section, we dilute randomly our A^-body sample in two kind 
of subsamples (see TableO: 

(Dl) 125 realizations of dark matter catalogs involving 50000 
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: 2. Same as Table[T] but for the dilute dark matter sub-samples, 


Dl andD2. 










Samp. Content 


A'rca 


Size 






Measured f! from: 


Comment 












_u/__/_ 


/ 2 

























(Dl) 50000, real space 


125 


L = 200 


128 


50000 


0.163 


0.213 


0.278 


( all ), no smoothing 












0.004 


0.006 


0.008 


error from dispersion 












0.237 


0.252 


0.268 


( all ), smoothed 












0.007 


0.008 


0.009 


error from dispersion 












0.283 


0.300 


0.319 


( 1.5(7 isocontour ), no smoothing 












0.011 


0.012 


0.013 


error from dispersion 












0.277 


0.288 


0.299 


( 1.5(T isocontour ), smoothed 












0.012 


0.013 


0.013 


error from dispersion 


50000, redshift space 


125 


R = 70 


128 


8980* 


0.062 


0.248 


3.841 


( all }, no smoothing 












0.025 


0.027 


16.11 


error from dispersion 












0.100 


0.184 


0.378 


( all ), smoothed 












0.027 


0.022 


0.118 


error from dispersion 












0.227 


0.264 


0.309 


{ 1.5(T isocontours ), no smoothing 












0.030 


0.034 


0.050 


error from dispersion 












0.209 


0.237 


0.269 


( 1.5(T isocontours ), smoothed 












0.026 


0.028 


0.039 


error from dispersion 


(D2) 11000, real space 


125 


L = 200 


128 


11000 


0.192 


0.250 


0.326 


( all ), no smoothing 












0.011 


0.017 


0.026 


error from dispersion 












0.216 


0.255 


0.300 


( all ), smoothed 












0.015 


0.020 


0.026 


error from dispersion 












0.251 


0.292 


0.340 


{ 1.5(T isocontour ), no smoothing 












0.023 


0.027 


0.033 


error from dispersion 












0.249 


0.281 


0.318 


( 1.5(T isocontour ), smoothed 












0.024 


0.028 


0.033 


error from dispersion 


1 1000, redshift space 


125 


i? = 70 


128 


1976* 


0.074 


0.226 


1.323 


( all ), no smoothing 












0.028 


0.035 


1.687 


error from dispersion 












0.096 


0.197 


0.492 


( all ), smoothed 












0.031 


0.034 


0.238 


error from dispersion 












0.190 


0.253 


0.347 


( 1.5(T isocontours ), no smoothing 












0.040 


0.048 


0.099 


error from dispersion 












0.182 


0.233 


0.308 


( 1.5(7 isocontours ), smoothed 












0.040 


0.044 


0.089 


error from dispersion 



points, similarly as for the mock with mass thresholding Mth ~ 
5.10" Mq used in §1331 

(D2) 125 realizations of dark matter catalog involving 11000 
points, similarly as for the mock catalogs with mass thresholding 
Mth = 4.10^^ A/0 used in § 1441 

We perform exactly the same measurements in these dilute catalogs 
than in sample (i) (in real space, except that we use 128'' grids 
for interpolating the fields and we have 125 different realizations) 
and in samples (iv) (in redshift space). In the latter case, note that 
the observer position in the simulation cube is different for each 
realization, chosen exactly to be on a regular patter covering the 
full simulation volume as in samples (iv). 

The results are summarized in Table |2l and illustrated by 
Figs. |7l and |8l We first discuss real space measurements using the 
second order moments of the full PDF and compare the values ob- 
tained for n to those from the undiluted sample (i), given in Table[Tl 
Note that this sample uses a 512'^ grid for the interpolation instead 
of the 128'^ one for samples Dl and D2. Using a 128'' resolution 
grid increases the measured value of Q in first line of Table[T]from 
0.157 to 0.178. Taking that fact into account, we notice that, due to 
the dilute nature of the samples, our interpolation procedure uses 
an adaptive kernel of larger size: as shown by Fig.|71 this makes the 
fields more linear, closer to the Gaussian limit, and decreases the 
level of effective bias brought by the "propeller" shape of the bivari- 
ate PDF. As a result, prior to additional smoothing with a Gaussian 



window of size = 10 Mpc, the measured Q is larger for the 
dilute samples than for the full one, and the convergence with linear 
theory prediction improves with level of dilution. The mean inter- 
particle distance in the sparser sample is of the order of A = 9 /i^^ 
Mpc, slightly lower than the size of the post-processing Gaussian 
smoothing kernel. One thus expects rough agreement between full 
sample and samples Dl and D2 after smoothing with such a win- 
dow, which is approximately the case. Furthermore, the measured 
value of f2 in D2 is not very sensitive to whether additional smooth- 
ing is performed or not, since A ~ ^. 

Still examining real space measurements, we now consider the 
values of measured from second order moments of the PDF, but 
using only the region of best likelihood, Eq. J52b . similarly as in 
previous paragraph. As expected, the effective bias due to nonlinear 
contributions (or propeller effects) is tremendously reduced, and 
one recovers a value of O compatible with the true value, given 
the errorbars. These latter, which estimate pure shot noise, are of 
the order of 4 and 10% for Dl and D2 respectively. As for the full 
sample, this result stands also when additional smoothing is applied 
to the data, although one can notice a slight bias to lower values of 

n. 

We now turn to redshift space measurements and first consider 
the measured value of Q using the second order moments of the full 
PDF. The results obtained in previous section still hold (compare 
Table [il to Table |2l(, except that the measurements performed be- 
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RAMSES 50000, Real space, Lj = 1.6h"' Mpc 



RAMSES 50000, Real space, L=10h"' Mpc 
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RAMSES HOOO, Real space, L=10h"' Mpc 
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Figure 7. Similarly as in Fig. [3] but for the 125 realizations of diluted samples Dl (upper panels) and D2 (lower panels), the left panels con'espond to the 
raw interpolated fields on a 128"^ grid, while for the right panels, a smoothing with a Gaussian window of size = 10 h^^ Mpc was performed prior to the 
measurements. Since the samples are now more diluted, the adaptive interpolating kernel is of much larger size, which decreases considerably nonlinear tails: 
the propeller shape of the PDF is less apparent but it still affects the measurements: the slope of the dotted line is larger than the linear expectation given by 
the thick solid line. 



fore additional Gaussian smoothing give a lower value of fi, 0.248 
for Dl and 0.226 for D2: the low bias effect on SI gets worse with 
dilution, since the larger size of the adaptive kernel tends to reduce 
the effect of fingers of God: as noticed previously, finger-of-God ef- 
fects help to reduce the asymmetry brought by the "propeller", the 
main source of the low bias on Q,. Additional Gaussian smoothing 
with a window of size ^ = 10 h^^ Mpc improves the convergence 
between Dl, D2 and the full sample, as expected, but induces a 
highly underestimated value of SI ~ 0.2, because of the propeller 
effects which are then prominent. These arguments are supported 
by examination of Fig. [8] Note that the asymmetry of the tails of 
the bivariate PDF around the major axis of the maximum likeli- 
hood region, already observed and explained in Fig. (6] is now more 
pronounced, at least from the visual point of view. 

Measurements are improved while selecting the region of best 
likelihood, but not as well as in the real space case or in the red- 
shift space case with full sampling: at best, Q, is underestimated 
by about 12 percent. Additional smoothing or passing from Dl to 
D2 expectingly increases the bias. This underestimate comes again 
from the fact that the adaptive interpolating kernel is now of much 
larger extension than for the full particles sample, which induces 



biases comparable to what was observed for the full samples with 
additional smoothing (ij |4.3.3l l. 

Note that it is possible to reduce the effective bias observed 
on S7 by narrowing the region of best likelihood at the price of 
an increase of the errors. For instance, taking a 38 percent confi- 
dence region enclosed by a ~ lo isocontour gives SI = 0.285 
and Q, — 0.273 for Dl and D2, respectively. Our procedure for 
measuring Q can thus certainly be improved with a more sophis- 
ticated treatment of the region of best likelihood. For instance, a 
way to extract in an unbiased way the parameter Q from the data 
could consist in measuring the local slope of the skeleton of the 
surface representing the bivariate PDF (see Novikov, Colombi & 
Dore, 2006) after appropriate (adaptive) smoothing of the velocity- 
gravity scatter plot. This is left for future work; in what follows, we 
shall still use for the sake of simplicity our l.Scr likelihood contour 
technique, while staying aware of the bias brought by dilution. 



4.4 Measurements on simple mock "galaxy" catalogs 

To estimate in a sufficiently realistic way how biasing affects the 
results on the velocity-gravity relation, we extracted from the sim- 
ulation four "galaxy" catalogs, corresponding to two methods of 
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Figure 8. Similarly as in Fig.|6] but for the diluted samples Dl (upper panels) and D2 (lower panels). 



treating dark matter halos. In the first method, we consider each 
dark matter halo as a galaxy. In the second method, we consider 
each substructure present in dark matter halos as a galaxy. The de- 
tails are given in § 14.4.11 For each kind of "galaxy catalogs", two 
level of biasing are considered, corresponding to small and large 
threshold on the dark matter (sub)-halos masses, M > 5.10" Mq 
and Af > 4.10^^ Mq respectively. In all cases, each galaxy is 
given the same weight while density and velocity are interpolated 
as explained in § 14.11 This rather extreme way of treating the bias 
is far from being optimal: appropriate weighting could be given to 
each object in order to correct at least partly for the effects of the 
bias, as further discussed in the conclusions of this paper. 

The quantitative analysis of the velocity-gravity relation is 
performed in §|4A2] The measurements are interpreted in the light 
of the previous paragraphs analyses, using as a guideline the large 
scale bias obtained from the measurement of the power-spectrum 
of the "galaxy" distribution. 



Colombi, 2004)0 AdaptaHOP builds an ensemble of trees. Each 
tree corresponds to a halo which is a connected ensemble of par- 
ticles with SPH density p > 81. The branches of the trees are 
composite structures of which the connectivity is controlled by the 
saddle points in the particle distribution. The leaves of the trees 
correspond to the smallest possible substructures one can find in 
the simulation. From this ensemble of trees, we extract two kinds 
of catalogs, one where each galaxy is identified to a tree, the other 
one where each galaxy is identified to the leaves of the trees (if a 
tree has only one leaf, it means that the halo is its own single sub- 
structure). Note that velocities of these galaxies are computed as the 
average velocity of all the particles belonging to the corresponding 
(sub-)structure. 

Additional mass thresholding is used to control the number of 
"galaxies" in the catalogs,. In final, 4 mock catalogs are obtained 
(see also TableO: 

(v) A "halos" catalog, involving 43482 dark matter halos with 
masses larger than Afth = 5.10^^ A/q; 

(vi) A "massive halos" catalog, involving 11934 dark matter ha- 
los with masses larger than A/th = 4.10^^ Mq; 



4.4.1 The mock catalogs 

To extract halos and substructures from the simulation, we use 
the publically available software adaptaHOP (Aubert, Pichon & 



° The parameters used in adaptaHOP are the same as in Aubert et al. 
(2004, Appendix B), namely A^sph = 64, A^hop = 16, Pth = 81 
and /poiason = 4. 
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Table 3. Same as Table[T] but for the "galaxy" catalogs, (v)-(viii). In the second column, it is specified if they coirespond to a halo ("HAL.") or a substructure 
catalog ("SUBS") and what mass threshold was used to select the dark matter (sub-)structures. The values of O displayed in columns seven to nine are obtained 
from the measured /3 assuming linear theory predictions with no bias. As a guidance to understand the results, we show in the second column of the table the 
measured value of the bias b at the simulation box size scale (using the power-spectrum of the interpolated density), as well as f2(,fl = fl/b^^^, the expected 
value of f2 from the measured /3 if unity bias is assumed. 



Samp. 


Content 




Size 








Measured Q from: 


Comment 


(v) 


HAL. > 5cllM0, real space 


1 


L = 


200 


128 


43482 


0.304 


0.375 


0.464 


all, no smoothing 




b = 0.80, floH = 0.45 












0.330 


0.379 


0.435 


all, smoothed 
















0.318 


0.368 


0.426 


1.5(T isocontour, no smoothing 
















0.329 


0.370 


0.415 


1.5(T isocontour, smoothed 




HAL. > 5ellM0, red. space 


125 


R = 


70 


128 


7874* 


0.208 


0.313 


0.487 


( all ) , no smoothing 




b = 0.80, Qoff = 0.45 












0.033 


0.035 


0.089 


error from dispersion 
















0.259 


0.312 


0.378 


( all ) , smoothed 
















0.035 


0.039 


0.055 


eiTor from dispersion 
















0.311 


0.360 


0.419 


( 1.5(T isocontours ), no smoothing 
















0.054 


0.063 


0.082 


eiTor from dispersion 
















0.310 


0.345 


0.385 


{ i.Scr isocontours ), smoothed 
















0.058 


0.064 


0.078 


error from dispersion 


(vi) 


HAL. > 4el2M0, real space 


1 


L = 


200 


128 


11934 


0.206 


0.254 


0.313 


all, no smoothing 




b = 0.95, QaS = 0.33 












0.218 


0.254 


0.296 


all, smoothed 
















0.225 


0.259 


0.297 


i.Scr isocontour, no smoothing 
















0.231 


0.258 


0.288 


i.Scr isocontour, smoothed 




HAL. > 4el2MQ, red. space 


125 


R = 


70 


128 


2143* 


0.155 


0.220 


0.320 


( all ) , no smoothing 




b = 0.95, fieff = 0.33 












0.029 


0.030 


0.055 


eiTor from dispersion 
















0.178 


0.221 


0.277 


( all ) , smoothed 
















0.030 


0.032 


0.045 


error from dispersion 
















0.215 


0.255 


0.305 


( 1.5(T isocontours ), no smoothing 
















0.044 


0.049 


0.065 


error from dispersion 
















0.215 


0.245 


0.282 


{ i.Scr isocontours ), smoothed 
















0.045 


0.049 


0.061 


error from dispersion 


(vii) 


SUBS. > 5ellA/Q, real space 


1 


L = 


200 


128 


55044 


0.215 


0.253 


0.298 


all, no smoothing 




b = 0.97, QcB = 0.32 












0.248 


0.265 


0.283 


all, smoothed 
















0.260 


0.281 


0.303 


1.5cr isocontour, no smoothing 
















0.264 


0.277 


0.292 


1.5cr isocontour, smoothed 




SUBS. > 5ellA/Q, red. space 


125 


R = 


70 


128 


9886* 


0.101 


0.228 


0.607 


( all ) , no smoothing 




b = 0.97, Oeff = 0.32 












0.024 


0.020 


0.223 


eiTor from dispersion 
















0.144 


0.198 


0.277 


( all ) , smoothed 
















0.023 


0.019 


0.041 


error from dispersion 
















0.234 


0.264 


0.299 


( 1.5fT isocontours ), no smoothing 
















0.030 


0.035 


0.045 


error from dispersion 
















0.224 


0.244 


0.267 


( 1.5cr isocontours ), smoothed 
















0.028 


0.031 


0.038 


error from dispersion 


(viii) 


SUBS. > 4el2MQ, real space 


1 


L = 


200 


128 


11221 


0.189 


0.226 


0.269 


all, no smoothing 




b = 1.035, Qoff = 0.28 












0.204 


0.228 


0.254 


all, smoothed 
















0.215 


0.239 


0.266 


1.5cr isocontour, no smoothing 
















0.218 


0.236 


0.255 


1.5cr isocontour, smoothed 




SUBS. > 4el2M0, red. space 


125 


R = 


70 


128 


2015* 


0.102 


0.176 


0.327 


( all ) , no smoothing 




b = 1.035, Qcff = 0.28 
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(vii) A "sub-structures" catalog, involving 55044 dark matter 
sub-halos with masses larger than Afth ~ 5.10^^ Mq; 

(viii) A "massive sub-structures" catalog, involving 11221 dark 
matter sub-halos with masses larger than Mth ~ 4.10^^ Mq. 

From each of these catalogs, we compute the velocity and gravity 
field on a 128^ grid as explained in §|4T| in real and redshift space. 
In the latter case, we use 125 different observers and perform the 



measurements in a sphere of radius 70 Mpc Mpc centered on 
the observer exactly as was done to generate the pure dark matter 
realizations (iv) in 14.31 

The philosophy underlying to the making of these catalogs 
relies on the facts that: 

(a) dark matter substructures are reasonable tracers of galaxies 
(e.g. Springel et al. 2001; Weinberg et al. 2006). In that sense, one 
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can consider the sub-structure catalogs as the most realistic. This 
is known to be true only to a limited extent. Indeed, sub-halos tend 
to be tidally disrupted while they spiral in their host halo, therefore 
one expects less sub-structures than galaxies in the core of dark 
matter halos (e.g. Diemand et al. 2004; Nagai & Kravtsov 2005). 
This is all the more true since we apply a mass thresholding to 
control the number of objects in our catalogs. However, the effect 
of such a depletion should be noticeable only at the smallest scales 
and should not affect significantly the results of the analyses at the 
level of accuracy reached in this paper. 

(b) dark matter halos are representative of the galaxy distribu- 
tion in terms of an ensemble of structures composed of clusters, 
groups of galaxies and field galaxies. Therefore, if one considers 
the galaxy distribution from a slightly different perspective, these 
structures can be used as well to study the velocity-gravity relation, 
with the appropriate weighting. The potential advantage of such an 
approach is to reduce considerably the fingers of God effects dis- 
cussed previously, since all the galaxies belonging to one cluster or 
group of galaxies are collapsed to a single point. Note that we over- 
simplify the analyses here by purposely giving the same weight to 
all the dark matter halos: this rather extreme procedure is expected 
to introduce quite significant (anti-)biasing effects on the gravity 
field determination. 



4.4.2 Quantitative measurements and the ejfect of the bias 

The results of our analyses are summarized in Table [3] and illus- 
trated by Figs. |9] and [To] As shown in 14.31 due to nonlinear ef- 
fects, measurements based on the direct estimation of the moments 
of the joint probability of velocity and gravity tend to bias the es- 
timated linear 13 to lower values, both in real and redshift space. 
To reduce this systematic effect, we propose to consider only the 
region of best likelihood. Strictly speaking, this approach makes 
sense only if the bias between the galaxy distribution and the dark 
matter distribution is either inexistent or linear. The concept of a 
non local, scale dependent bias as experienced here complicates 
considerably the analyses. However, the aim is to capture the lin- 
ear contribution of the dynamics, which in terms of our volume- 
weighted measurements is dominated by underdense regions. In 
these regions, the size of the adaptive kernel is large: inside the 
contour of best likelihood of the joint PDF, the measurements are 
expected to be dominated by the large scale bias, fe(fc), with small 
wavenumber k. In that regime, the bias is in general roughly linear 
or close to linear. The need to actually consider the region of best 
likelihood of the PDF to perform the measurements is therefore 
furthermore justified. As a result, in what follows, only results of 
Table |3] obtained from selecting the region inside the 1.5cr contour 
of the PDF are considered. Since additional smoothing tends to in- 
troduce additional bias in the redshift space measurements, it will 
not be discussed either, although the results are shown in Table [3] 
for completeness. 

We first discuss real space measurements and try to under- 
stand the effects of the bias. Sample (v) is the halo catalog with 
low mass threshold. This catalog is expected to present signifi- 
cant antibias. Indeed, if small halos, which mainly lie in moder- 
ately dense and underdense regions, trace rather well the underly- 
ing dark matter distribution, the largest ones, which are collapsed 
to a single point, induce a significant underestimate of the over- 
all strength of the dark matter gravity field, resulting in the antib- 
ias effect. As shown in Table [3] the value of the bias measured 
at the simulation box size (using the Fourier modes at the largest 
scales) gives b{k = 2ti/L) = 0.8. If one uses, as argued above. 



this value of 6 as a reference, the effective value of Q, obtained 
from linear theory, ignoring the bias, is larger than the true one, 
f^eff = Q./h'^^^ = 0.45. This is to be compared to the value of 0.37 
obtained from ^y{g^/{v^■ Although this latter result is signifi- 
cantly lower than the "expected" value, two facts have to be taken 
into account: 

(a) As an effect of non local, scale dependent biasing, the cor- 
relation between gravity and velocity field is not as tight as for the 
dark matter diluted sample, DI, albeit this latter contains the same 
number of objects. Indeed, one measures SI = 0.31 and 0.42 from 
{g^)/{vg) and {vg) / {v'^) respectively in sample (v), while Q, was 
ranging between 0.28 and 0.32 for Dl. This loss of tightness in the 
velocity-gravity relation obviously affects in a non trivial way the 
region of best likelihood. 

(b) Our approach for estimating the effects of the bias is rather 
crude and does not give account of the additional subtleties related 
to our adaptive smoothing procedure. 

Taking these two points into account, the measurements are rather 
consistent with the expectations, at least at the qualitative level. 

In the halo catalog (vi), small mass "galaxies" are removed, 
so the underdense regions are not populated anymore, leaving us 
with a set of more clustered objects than in sample (v), since the 
higher is their mass, the more significantly the halos are clustered. 
The consequence is that the antibias found in catalog (v) is re- 
duced to b{k = 2ti/L) — 0.95, resulting in an effective ex- 
pected value of r^cff ~ 0.33, to be compared with the value of 
0.26 obtained from the measurement of •\/ {g^}/{v^). Again the 
arguments (a) and (b) developed just above apply, leading us to 
conclude that this result is consistent with expectations, and sim- 
ilarly for the substructures catalogs: sample (vii) is the closest to 
the underlying dark matter distribution, with a very slight antib- 
ias of b{k — 2n/L) — 0.97 and il^s = 0.32 to be compared 
to 0.28 from •\/ (g^) / {v^)\ sample (viii) presents a small positive 
large scale bias, with b{k — 2n/L) — 1.035 and Q^s ~ 0.28, to 
be compared to 0.24 from \/ 

We can conclude here that we understand what are the var- 
ious effects influencing the measurements in the mock catalogs. 
The general trend is that the "true" value of parameter (3 (i.e. the 
very large scale one) seems to be underestimated by about 10 to 
20 percent. We give here quotes, because such a value is not well 
defined. Indeed, the interpretation of the measurements is compli- 
cated by adaptive filtering, combined with a non local, scale de- 
pendent bias. But as discussed above, we really want to capture the 
value of /3 in the linear regime limit, where the bias is expected 
to be roughly linear, scale independent and in fact rather close to 
unity, as observed in our "galaxy" catalogs. Note interestingly that 
the conditional moment (vg) / (u^) gives a rather good estimate of 
the expected value of (3. Without going too far exploiting such a 
property, which may be specific to our catalogs, it is clear that the 
slopes given by the two conditional moments could be used as prox- 
ies to estimate a range of possible values for (3, or in other words, 
errorbars. Indeed, the difference between these slopes is related to 
the width of the likelihood region, i.e. it estimates the tightness 
of the correlation between gravity and velocity as measured in the 
catalog. However, note that this difference does not properly give 
account of systematic effects discussed at length in this paper, as 
well as cosmic variance effects, although it is indirectly related to 
them. 

We now turn to redshift space measurements. What matters is 
that they should be self-consistent with the real space ones. This is 
nearly the case, except again that they tend to underestimate real 
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Figure 9. The velocity-gravity relation in real space as measured in our mock "galaxy" catalogs (v)-(viii). The left column of panels corresponds to the direct 
SPH alike interpolation on a 128'^ grid, while an additional Gaussian smoothing with a window of size 10 h^^ Mpc was performed in the right column of 
panels. The informations on the catalogs are given as titles on each panel, namely, from top to bottom, halo catalog with mass thresholding M^^ = S.IO^^Mq, 
sub-structure catalog with Afth = 5.10^^A/q, halo catalog with M^^ = 4.10^^A/q and sub-structure catalog with M^^ = 4.1O^^Af0. The prediction 
from linear theory in absence of bias is given as the solid line while the slope obtained directly from the ratio (g| ) / (v| ) is represented as a dotted line. 



space measurements by up to ten percent, an effect of sparse sam- 
pling discussed in ij |4.3.4l with no additional surprise. Overall, the 
difference between the measured value of Q in redshift space and 
the "expected', effective one, Slcff , is about 25 percent, with a sys- 
tematic bias to lower values, as always. Obviously our mock cata- 



logs are very peculiar in a sense that they represent quite an extreme 
case in the range of possibilities. And still, without any assumption 
on the bias, or in other words, assuming that b is of the order of 
unity, we can determine Q, with an accuracy of the order of 30 per- 
cents without significantly strong prior. Furthermore, the quality of 
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Figure 10. Same as in Fig. |9] but in redshift space. The scatter plots are obtained from the average of 125 observer positions disposed on a regular pattern 
spanning the simulation volume. As expected, finger of God effects are more prominent on the sub-structures catalogs than on the halos catalogs. 



this estimate can certainly be improved with appropriate weighting 
of the galaxies. However, as discussed in the conclusion that fol- 
lows, the sources of systematic errors related to the instrumental 
noise or to the method used to provide velocities were not consid- 
ered here. 



5 CONCLUSIONS 

In cosmic density-velocity comparisons, since the density field is 
given solely in redshift space, redshift-space effects have to be ac- 
counted for. One method is to perform the analysis directly in red- 
shift space, as proposed by ND (Nusser & Davis 1994). Another 
one is to reconstruct the real-space density field, but this method is 
complicated by the fact that redshift-space corrections require as- 
suming a value for p. As /? is not known a priori, one is forced to 
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reconstruct the density field for a range of values of which makes 
the real-space comparison involved, despite its apparent simplicity. 
In the present paper we have proposed a method that works in red- 
shift space, but is simpler than the formalism of ND, which is based 
on a spherical harmonics expansion. It is essentially as simple as in 
real-space; the only difference is that here a directly estimated pa- 
rameter is not just but its simple function. 

This method relies on a tight correlation between cosmic grav- 
ity and cosmic velocity in redshift space. We have derived the 
velocity-gravity relation in redshift space and its scatter analyti- 
cally, under the assumptions of the linear regime, distant observer 
limit (DOL) and CMB rest frame for redshift measurementsQxhis 
relation has turned out to be a simple modification of the corre- 
sponding one in real space and in the linear regime. Then we have 
performed a dark matter A'^-body experiment to test the velocity- 
gravity relation both in real space and in redshift space. We also 
extracted mock galaxy catalogs from this data set, with rather ex- 
treme prescription for selecting "galaxies" from the dark matter dis- 
tribution in order to analyze the effects of non-trivial biasing on the 
velocity-gravity relation. The main results are the following: 

• To perform the measurements, we propose to interpolate the 
gravity and the velocity fields on a regular grid using a new adap- 
tive smoothing procedure described in § 14.11 This method is close 
to SPH interpolation, but without the defect of having regions of 
space where the velocity field is undefined. Although our algorithm 
is expected to be less optimal than those based on Voronoi or Delau- 
nay interpolation (e.g., Bernardeau et al. 1997; Pelupessy, Schaap 
& van de Weygaert 2003), it has the advantage of being simpler to 
implement and much faster, and can be applied to very large sets 
such as the 512^ sample analyzed in this paper. 

• In real space, simple measurements using moments of the full 
joint probability distribution function (PDF) bias Q. to low values. 
This is due to the fact that these moments are contaminated by the 
tails of the PDF. These latter are influenced by nonlinear dynamics, 
which induces a global 'propeller' shape for the PDF. 

• In redshift space, fingers of God (FOG) induce additional 
"anti-diagonal" tails to the PDF. This effect seems to somewhat 
compensate the bias due to the propeller shape. However this result 
could be coincidental and is too sensitive to the level of smoothing 
applied to the data. Note that in actual peculiar velocity catalogs, 
clusters of galaxies are in general collapsed into one point to im- 
prove the signal to noise ratio of the distance estimates. Although 
FOG are thus less of an issue in real data, rich filaments can still 
contribute to an effect similar to FOG as discussed in i; 14. 3. II 

• Selecting about a l.So" (68% confidence) region around the 
maximum of the joint PDF of gravity and velocity fields makes the 
data to match the linear theory predictions very well, both in real 
and in redshift space. This result stands even in the highly nonlinear 
regime (no additional smoothing except the adaptive interpolation 
on the grid). It is easily explained by the fact that our measurements 
are volume weighted, giving most of the statistical weight to under- 
dense regions. Therefore, it is sufficient to exclude only the tails of 
the distribution to obtain unbiased estimate of Q,. 

• Additional smoothing does not improve the agreement with 
linear theory predictions: first, it is expected to increase deviations 
from the DOL limit; second, it makes the situation worse in the 
region of best likelihood, where it mixes linear contributions with 



nonlinear ones and biases 13 to lower values, an effect which is min- 
imized as much as possible with our adaptive smoothing procedure. 

• The measurements obtained in our mock catalogs provide a 
value of 13 in agreement with intuition, namely that large scale bias 
dominates if the region of best likelihood is selected. This is in- 
herent to our volume weighted approach which gives more weight 
to underdense regions (see also Berlind, Narayanan & Weinberg 
2001). However, due to the dilute nature of real catalogs, we ex- 
pect an effective bias on (3 to lower values. This bias was found 
to be of the order of 10-20 % for the mock samples considered 
in this paper. In addition, the velocity-gravity relation looses some 
tightness due to effects of non-local biasing between luminous and 
dark matter distributions. Nevertheless, our analyses show that di- 
rect measurement of f2 from real catalogs assuming no bias should 
give an answer accurate to about 30 % at worse. Moreover, this er- 
ror can probably be significantly improved to values as low as 10 
%, with the appropriate weighting of the data (purposely not per- 
formed here), as discussed below. 

There are several points we did not address in our numerical 
analyses, and we leave them for future work: 

- Noise: through the analysis of dilute samples and mock cata- 
logs, we addressed to some extent the effects of shot noise for real- 
istic modern catalogs. However, we used the same catalog for ve- 
locity and gravity fields, implying that the same adaptive kernel can 
be used for both fields, which is not realistic|f| Furthermore we did 
not take into account the fact that there is in practice a significant 
relative error on galaxy distance estimates, which can be as large as 
20 percent. Thus the velocity-gravity relation is certainly not ex- 
pected to be as tight in real catalogs as in the samples considered in 
this paper. Note interestingly, though, that since the measurements 
are performed in redshift space, issues related to Malmquist bias 
(e.g., Strauss & Willick 1995) should be irrelevant, as long as the 
distance estimates are unbiased (just noisy). 

- Incompleteness: Another closely related issue is that the cata- 
logs used to estimate the gravity field are incomplete. First, edge 
effects are expected to be significant because gravity force is of 
long range. The galaxy catalog used to estimate the gravity field 
has to be significantly deeper than the region used to perform the 
velocity-gravity comparison. We did not estimate in this paper how 
deep it has to be, not to mention the problem of obscuring by our 
own Galaxy. Second, the apparent density of galaxies decreases 
with increasing distance from the observer. In order to estimate the 
interpolated density field, one can, under the assumption of no bias 
between the galaxy distribution and the mass distribution, weight 
the galaxies by the inverse of the selection function, (f). In the linear 
regime, in the absence of bias and for redshifts measured relative to 
the CMB rest frame, the redshift space density contrast measured 
that way reads (Kaiser 1987) 

5,^5-12+—-^]--—. (54) 
\ amr / r or 

Under the DOL, the second term on the right hand side of this equa- 
tion drops, so the effect of the selection function disappears. How- 
ever, since remote parts of the catalogs are sparser, they might be 
given lesser statistical weight, hence augmenting the relative con- 
tribution of this deviation from the distant observer limit. 

- Better handling of the bias: our analyses in the simple mock 
catalogs show that biasing affects the velocity-gravity relation in 



Using CMB rest frame is crucial to avoid the so-called rocket effect on 
the redshift space gravitational acceleration, as discussed in §[5] 



° If the catalogs are different, a common smoothing kernel must be deter- 
mine to find a compromise between the velocity and the gravity samples. 
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a non trivial way. However, we used the most naive prescription 
to compute the gravity field, giving equal weight to all galaxies. 
Clearly, a better description would be to assign to each galaxy a 
weight proportional to its supposed host halo mass, for example 
assuming constant mass to light ratio, M/L, or a more sophisti- 
cated weighting using a function M/L = f[L) derived from ob- 
servations or obtained from theoretical models of galaxy formation. 
That would allow one to correct to a large extent for the effects of 
biasing, at the cost of additional priors. Still, unless the galaxy cat- 
alog used to compute the gravity field is very deep, i.e. includes 
faint galaxies, even giving the proper mass to light ratio to each 
galaxy does not correct for the fact that underdense regions can be 
artificially underpopulated, especially far away from the observer. 
To tackle with that, one could add a background population of low 
mass objects that would give account of the missing mass in the cat- 
alog and that would act as a shielding effect on the gravity field (see, 
e.g., Phelps et al. 2006). Again, some strong priors are required to 
deal with such a background, in particular on its clustering proper- 
ties and how it correlates with the population already present in the 
real catalog. 

Note finally that while our approach presents the advantage 
of simplicity compared to ND (and most other methods), we do 
not expect it to be as accurate as of ND, since it relies on the dis- 
tant observer approximation and uses in its current form adaptive 
smoothing, which complicates the interpretation of the results. (On 
the other hand, we clearly showed that the region of best likelihood 
of the joint velocity-gravity distribution is dominated by the lin- 
ear regime prediction, even in the highly nonlinear regime, which 
presents a noticeable advantage.) Similarly, our method is not ex- 
pected to perform as well as sophisticated reconstructions of La- 
grangian nature, which try to minimize the Euler-Lagrange action 
(e.g., Peebles 1989; Shaya, Peebles & Tully 1995; Nusser & Bran- 
chini 2000; Phelps 2002; Phelps et al. 2006) or to solve optimal 
assignment problem (e.g.. Croft &. Gaztaiiaga 1997; Frisch et al. 
2002; Mohayaee et al. 2003; Mohayaee & Tully 2005), although 
this remains to be verified. 
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APPENDIX A: SPHERICAL TOP-HAT IN REDSHIFT 
SPACE 

The general result in real space is 



9 



(Al) 



where 5 and 9 are the real-space density contrast and velocity diver- 
gence, respectively. In redshift space, in the distant observer limit, 
the general result is 
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In the linear regime, C = 2/3 and 9 = 5. Hence, 
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We see thus that the spherical top-hat model does not give the sta- 
tistical result obtained from linear theory. The reason for this is the 
following. In the spherical model, the velocity-gravity relation is 
fully deterministic, but we could also derive it calculating the rele- 
vant statistical averages. While in the linear regime, in general 

5l«'(k)=fe(l 



[where /i = k ■ s/{ks)], in the spherical top-hat we have 



(A7) 



= b{l+ (3/3)5, 



(A8) 



(A9) 



that immediately yields 

Si%{\^)=b{l + f3{^.'))5{k), 

since (/i^) = 1/3. In other words, the expression for the redshift 
space density contrast in the top-hat model is an average of gener- 
ally valid expression iAH over possible orientations of the vector 
k. In the case of the estimator (wsijs) / (vs), the average over angles 
yields an expression proportional to 



{^,^l + (3^l^)) 



i + p 



In the top-hat model, this is modified to 



{^^^l + ^3{^.'})) 



1 + 13 



(AlO) 



(All) 



Loosely speaking, the top-hat model performs a part of this average 
'too early'. So the top-hat model yields quantitatively different re- 
sults from the results of a rigorous statistical calculation. However, 
it is sufficient for qualitative purposes, needed here. In particular, it 
correctly predicts that the slope of the velocity-gravity relation in 
redshift space gets steeper. 

Turn-around in real space corresponds to extreme flatness of 
the pancake in redshift space. When 9 ~* 9t = 3/3^^, then g ^ 
and C tends to its maximal value, 2. Hence, 
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3/3-^ ( 



9t 3 



11. 



(A12) 



Here, 5t is related to 9t by the formula of Bernardeau (1992). In 
fact, Bernardeau et al. (1999) invented its more accurate modifica- 
tion: 



1 + 5- (l + e/a)° 



(A13) 



where a is slightly greater than 3/2. However, for our purposes 
Q = 3/2 should be accurate enough; moreover, it better describes 
the evolution of voids [for this value of a, 5(6' = —1.5) = — 1]. 
This formula works for mildly non-linear densities. 

For 9 > 3/3~^, the structure in redshift space is inverted: ve- 
locity, and so the ratio of gs to vg, changes sign. Then 



g = 139/3 
and 

93 _ 3 



1, 



C 



5 (3 



(A14) 



(A15) 



^^3 2 

For 9 = 9i = 6(3^^ we have g = 1, so the structure in redshift 
space is momentarily an (inverted) sphere, hence without further 
effort, C = 2/3. Therefore, then 

, 13^ 



V3 



-5.42. 



(A 16) 



For 9 ^ 3(3 ^ (highly nonlinear infall), the (inverted) structure in 
redshift space becomes very elongated (finger of God). Then g = 
(39/3 - 1 ~ (39/3 > 1 and e ~ {l/3)(39i. After some algebra, 
one can show that then 



C 



In {(39/3) 



{(39/3)^ ■ 

Moreover, then one cannot use formula l lA13t , but 5 
(Bilicki & Chodorowski, in preparation), hence 5/9 oc 
yields 
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(A 17) 

{9/3)' 
9. This 



(A 18) 



In other words, gs becomes very small compared to |u3| 
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